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In this note, which is of general stability theory interest, we discuss some of the key assertions 
usually stated in the context of the transition to turbulence problem. In particular, the two main 
points made here in the setting of the transition problem are (i) the crucial dependence of the 
stability results on whether the problem is considered on infinite or semi-infinite domain, and (ii) 
the energy conservation by the nonlinear terms of the Navier-Stokes equations. As an application, 
we demonstrate that the Couette flow, when analyzed in the mathematical setting most correctly 
reflecting the way the experiments are usually done, is spectrally unstable for finite Reynolds num- 
bers in apparent contradiction to the commonly accepted classical century-old results. Also, the 
interrelation of various stability notions, the effects of infinite dimensionality, the covariant nature 
of the transition phenomena and how non-normality of the linear operators and finite-amplitude 
instability fit into this picture are discussed as well. 
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I. INTRODUCTION 

When exploring physical phenomena, one naturally 
uses both experimental and theoretical tools to gain as 
complete understanding as possible. If understanding of 
the physics is relatively complete, then not fully rigor- 
ous mathematical description usually suffices. If, how- 
ever, the phenomenon is as complicated, as the transition 
to turbulence problem, where despite numerous experi- 
mental studies there is no firm grasp of physics, then 
the progress inevitably relies upon rigorous mathemati- 
cal tools. Since, all the equations we use in physics are 



phenomenological, then only a constructive symbiosis of 
experimental and rigorous theoretical tools allows us to 
understand the validity and the range of applicability of 
these equations. The Navier-Stokes equations (NSEs) for 
incompressible fluid, we discuss here, 
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are not an exception, since their analytical properties 
are still not fully understood. System (1) is supplied, 
as usual, with the appropriate boundary conditions at 
dfl and initial conditions at t = 0. 

The transition to turbulence problem has long been 
considered in the realms of hydrodynamic stability the- 
ory, the key question of which can be formulated as fol- 
lows: what happens to a given fluid flow (base state) 
under the influence of disturbances. If the flow is ro- 
bust under the influence of all possible disturbances, it is 
called (Lyapunov) stable and can be expected to be ob- 
served in Nature. If there are perturbations which start 
to grow, the flow is called unstable and thus is expected 
to break up. If there is some finite critical amplitude of 
disturbance, beyond which the flow is unstable, then it is 
called finite-amplitude unstable. One can also enrich this 
picture by quantifying not only amplitudes of the critical 
disturbances, but also their geometry in some appropri- 
ate infinite-dimensional phase space. Despite this seem- 
ingly straightforward view of stability, there is still no 
theory which would predict robustly the experimentally 
observed behavior, especially for canonical flows such as 
Couette, Poiseuille and Hagen-Poiseuille flows, where lin- 
ear spectral stability analysis is known to fail [1]. 

The current experimental evidence suggests that the 
transition from laminar to turbulent state has the char- 
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acter of a finite-amplitude instability, cf. the recent work 
of Mullin and Peixinho [2] and figure 7 in that reference. 
However, based on the natural limitations of any realis- 
tic experimental procedure, such as inability to introduce 
disturbances of all possible forms and to wait infinitely 
long to see if a particular disturbance will lead to tran- 
sition (the pipe is of finite length!), one has to keep in 
mind that the actual transition may turn out to be not a 
finite-amplitude instability [47], but an infinitesimal in- 
stability in a Lyapunov sense [48]. At the same time, 
due to the finite length of a pipe, one cannot say for sure 
that the disturbances, which seem to grow, will not decay 
eventually. 

Thus, it would be honest to say that the experimental 
suggestions cannot be considered as a firm evidence, and 
the gap between the natural deficiencies of experiments 
and the current theories can only be bridged by careful 
sorting out and making precise the theoretical results and 
by the subsequent punctilious interpretation. 

In this note we make one step in this direction by an- 
alyzing some of the key assumptions in the current theo- 
retical views. The first one is the standard consideration 
of disturbances as defined on domains which are infinite 
in the streamwise direction, — oo < x < +oo, which is 
usually justified by the translational invariance of the 
base plane-parallel state. As we will show in §11, the 
usage of a more relevant to experiments domains setup, 
> x < +00, leads to unexpected stability results. The 
second one lies at the foundation of the transient growth 
approach, namely the energy conservation by the nonlin- 
ear terms of the NSEs. Just to remind the reader, the 
arguments usually made in the context of the transition 
to turbulence problem are 

1. Nonlinear terms do not produce energy; 

2. Thus, in order to explain transition one has to focus 
on the linear terms; 

3. The linear terms produce energy only transiently 
and thus the transient growth, originating in the 
non-normality of the linear operator, is the key to 
explaining the transition, i.e. the transition is "es- 
sentially linear" [3] and non-normality is the neces- 
sary condition for subcritical transition [4, 5]. 

As we will show in §111, the energy conservation is not 
true in the context of the transition to turbulence prob- 
lems in general. In the rest of the paper, we discuss other 
important issues of the current theories, namely the in- 
terrelation of various stability notions, the effects of infi- 
nite dimensionality, the covariant nature of the transition 
phenomena and how non-normality of the linear opera- 
tors and finite-amplitude instability fit into this picture. 

II. DOMAIN TYPE EFFECTS 

As mentioned in the introduction, the stability analysis 
of the Couette and pipe flow is usually performed on do- 
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FIG. 1: Kovasznay flow. 

mains unbounded in both directions, — oo < x < +oo. 
If one recalls the way the experiments on the transi- 
tion are usually done, i.e. one introduces disturbances 
at the inlet location and observes how they evolve down- 
stream, then it becomes clear that the semi-infinite do- 
main, x G [0, +oo), is more relevant as a mathematical 
idealization (in reality the domains are finite, of course) . 
To illustrate the domain type effect on the stability re- 
sults, namely whether the domain is infinite or semi- 
infinite, we consider first the Kovasznay flow in §11 A, 
which is treatable analytically and thus allows us to make 
the key points in the most clear way, and then the Cou- 
ette flow in §11 B, which we study both numerically and 
analytically. 



A. Kovasznay flow 

Let us first analyze the inviscid version of the Kovasz- 
nay flow, i.e. the flow behind a periodic grid located at 
x = 0, as shown in figure 1. This choice of the flow is 
due to its direct relevance to the transition to turbulence 
problem both in terms of the plane-parallel base state 
and the flow domain type, as well as due to transparent 
analytical treatment. In the 2D stream- function formu- 
lation, when the velocity components are u = ip y and 
v = —ipxt its dynamics obeys 

A^ t + ip y A^ x - V*AV>j, = 0, (2) 

which after decomposition, tp = "J - + ip', into the basic 
state &(x, y) = y periodic for x > with period 1, and 
into the perturbation ip' produces the following linearized 
evolution equation for i/V: 

A$ + * V A< = 0. (3) 

Performing the standard eigenvalue analysis of (3), 
ip'(t,x,y) — > (f>(x,y)e xt , where in view of the solution 
periodicity in y-direction 

oo 

4>(x,y) = a n (x) cos 2irny + (3 n (x) sin 2nny, (4) 

n=0 
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FIG. 2: Semi- versus infinite flow domains. /+ is defined on 
x G [0, +00), /_ is defined on £ £ (— 00, 0], / is defined on 
x G (— oo, +00) 



we find that all the mode amplitudes a n (x) and f3 n (x) 
are decoupled, since the operator C = (A + d x ) A acting 
on <f)(x,y) is of even order in d y . Hence, 
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(5) 



that is each "amplitude" a n satisfies 

(A + d x )(dl - An 2 n 2 )a n {x) =0,n£ Z 4 



(6) 



For example, for n — let us contrast the solution of this 
eigenvalue problem on infinite and semi- infinite domains: 

• x G (— 00, +00): assuming that y(x) G L2, by ap- 
plying Fourier transform we get A = —ik, k G M, 
i.e. marginal stability. 

• x G [0,+oo) and \y(x)\ < 00: clearly, instability 
is present since there is an eigenfunction y ~ e _Ai2; , 
/i > such that the eigenvalue A = /x, which clearly 
leads to instability. 

Same type of analysis can be done in the viscous case as 
well. 

This counter-intuitive difference in the stability re- 
sults between the semi-infinite and infinite domains can 
be appreciated with the sketch in figure 2. Namely, 
if, for example, one restricts (eigen-) functions f(x) to 
be bounded for all x including infinities as motivated 
by the physics, then the space of functions defined on 
x G (—00, +00) is more restricted compared to the space 
of functions defined on x G [0, +00). Indeed, if one 
can construct a function /+ bounded on x G [0, +00) 
which also satisfies the Orr-Sommcrfcld (OS) equation 
(7a), then continuation of this function onto x G (—00, 0] 
may lead to an unbounded function /_, as dictated by 
the structure of the linear Orr-Sommcrfcld operator and 
as illustrated in figure 2. Note that this explains the 
sensitivity of the critical bifurcation (Reynolds) number 
to the properties of disturbances at the domain inlet, 



x = 0: while their amplitudes do not play a role in view 
of the linearity of the problem, gradient-like properties 
of the disturbances do! Indeed, varying these properties 
of disturbances at x = 0, which may be masked by their 
amplitude, effectively changes the boundary conditions 
at x — and thus the size of the function space. Since 
restricting the domain to x G [0, +00) enlarges the func- 
tion space, one can expect that the spectrum enlarges as 
well and may lead to instabilities. 

Finally, just to reiterate on the crucial distinction of 
the stability on infinite versus semi-infinite domain, we 
remind the reader that an instability on a semi-infinite 
domain implies a growth of some eigenmode, say function 
/+ defined only on a semi-infinite domain in figure 2. 
Then, if the domain is infinite, the unstable eigenmode 
/+ obtained for a semi-infinite domain is not defined in 
general (cannot be continued for negative x in the original 
function space, say the space of bounded function) and 
cannot grow in the part of the domain, where it is not 
defined. 



B. Couette flow 

One can expect that the analogous behavior should 
take place in the Couette and other "troublesome" flows. 
First recall that in the case of the Couette flow, Romanov 
[6] rigorously proved that the disturbance defined on x G 
(—00, +00) and decaying at x = ±00 can not lead to 
instability. 

Let us show that the Couette flow is unstable on a semi- 
infinite domain as opposed to the case considered on an 
infinite domain. In the latter case, when the disturbance 
is defined on— 00 < x < +00 and decays at infinities, 
it is known [6] that the upper bound of the real parts 
of the spectrum of the linear operator A is sup Re(A) < 
—c/Re with c > 0. If, on the other hand, analogously to 
the Kovasznay flow we consider the eigenfunctions </>(x, y) 
defined on x G [0, 00), then assuming the separated form 
(j>(x,y) = a(y)b(x) = a(y)e~^ x with fx > [49], we get 
the following eigenvalue problem for a(y): 
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yy ""uyyyt 
— 1,1: a = a y = 0, 



which is given in the case /1 = (which according to the 
Fourier analysis should be present if the disturbance does 
not decay at x = 00). The straightforward dispersion 
relation for this equation yields 

A = 47rV Re' 1 , n G N, 

i.e. there is the eigenvalue A = which gives marginal 
stability for any Re, as opposed to the result on the infi- 
nite domain. 

While the above is the rigorous result on the stability 
of the Couette flow when the disturbances do not decay 
at infinity, one can apply heuristic Syngc's method to the 
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general case, fx > 0: 



d 



— + /i - iie (A - ny) 
dy z 



^+M 2 )a(y) = 0, (7a) 
y = — 1, 1 : a = a y = 0. (7b) 



Multiplying by a*, integrating over t/e [—1,1], and using 
the notation if = |a^| 2 dy, we get 

Ee A(/i 2 / 2 - J 2 ) = if - 2 M 2 / 2 + M 4 / 2 

-i?e/i / yal&y-jj 2 \ ya 2 dy 
Jo J -i 

which in the limit [i — > +oo, physically corresponding to 
the localized disturbance at the inlet, yields 



(8) 



i.e. one should observe an instability. The above asymp- 
totics is, of course, valid only if the solution a(y) does 
not have a strong dependence on /i. While, as we will see 
shortly, a(y) docs depend on [i, (8) turns out to give the 
right insight. 

Rigorous dispersion relation for the Couette flow on 
semi-infinite domain can, in fact, be written down analyt- 
ically. If we rewrite (7) in the operator form, LiL 2 a(y), 
where L x = + - Re (A - fxy) and L 2 = + M 2 , 
then we can exploit this factorization of the original forth- 
order operator in deriving the dispersion relation. First 
note that the kernel of these operators without boundary 
conditions are 

ker L\ = span {Ai(. . .), Bi(. . .)} , 
ker L 2 = span {sin(. . .), cos(. . .)} . 



Thus, effectively we are solving 

Lfa(y) =aAi(...)+/?Bi(...), 



(9) 



where L 2 ° is the smaller operator with the boundary con- 
ditions (7b). Thus, in order to solve (9) one must have 

aAi(. ..)+/3Bi(. ..)_LkerL™*, 

where L9, * is the adjoint (not formally adjoint!) to L® . 
Simple usage of the definition of adjoint 



(v,L° 2 °w) = (Z%*v,w) = (L 2 v,w) 



(10) 



shows that L 2 °* coincides with L 2 , i.e. the operator with- 
out boundary conditions! Since we know the kernel of L 2 , 
then for equation (9) and thus (7) to have a solution, it 
is necessary and sufficient that the Gram determinant 
vanishes: 



(Ai| cos) (Ai| sin) 
(Bi|cos) (Bi|sin) 



= 0. 



In order to get some useful insights into the structure of 
a solution, let us analyze one of the entries of the above 
dispersion relation, e.g. 



(Ai| cos) 



Ai(— ay + (3) cosfiydy, 



where 
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As suggested by the asymptotics (8), it is tempting to 
consider the limit [i 3> 1 and thus to exploit the sta- 
tionary phase method, i.e. by considering cos \xy as fast 
oscillating function and Ai(— ay + (3) as a slowly varying 
function. However, it turns out that for the given range 
of y the periods of oscillation of both Airy and cosine 
functions are of the same order, which invalidates appli- 
cation of the stationary phase method. Another possible 
approach is to consider the Airy function with large ar- 
guments I — ay + /3| 3> 1, but then either A ~ O(l) leads 
to inconsistent asymptotics or |A| ^> 1 invalidates the as- 
sumption of a uniform large argument | — ay + [3\ ^ 1 
for all y. In any case, large values of [i indeed turn out 
to produce interesting structure of the solution, as illus- 
trated in figure 4 for [i; higher values of [i increase the 
number of oscillations in the "accordion" structure. 

In view of this fundamental difficulty to resolve this 
problem with the available rigorous analytical methods, 
we appealed to the numerical solution of (7) by expand- 
ing the solution into functions based on the complete set 
of the Chebyshev polynomials T n (y) = cos narccos y 



N 



(ii) 



i=0 



which guarantee the convergence faster than any power 
of N [7] . Here the basis functions are given by 



My) = T iiv) 



l+A 1+6 



T i+ i(y), 



i.e. they automatically satisfy the boundary conditions 
(7b). All the inner products associated with the Galerkin 
projection can be computed analytically, which is another 
advantage of this method. 

Since the goal here is not to study the complete bifur- 
cation picture for the Couette flow, but to demonstrate 
that it is spectrally unstable for some finite Reynolds 
number, say Re = 5000 and (i = 1. The resulting spec- 
trum is shown in figure 3(a) and the eigenfunction corre- 
sponding to the eigenvalue with the largest real part in 
figure 3(b). As figure 3(a) clearly suggests, the transi- 
tion in the Couette flow is of Hopf bifurcation type. For 
each fixed fi one can compute a critical Reynolds num- 
ber Re c , e.g. for /i = lwc get Re c ~ 75. As fi increases, 
i.e. the disturbance is localized around the inlet, the 
value of Re A max increases as well, as anticipated (8) and 
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(a) Spectrum (eigenvalues continue to 
the negative part of the real axis. 




05 1 , ,\i 

a(y) 

(b)Most unstable eigenfunction. 

FIG. 3: Spectral picture for the Couette flow at Re — 5000 
and \i = 1. The rightmost eigenvalue has real part ReA max ~ 
0.75988755 and the corresponding eigenfunction is shown in 
figure 3(b). 
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FIG. 4: The eigenfunction corresponding to the rightmost 
eigenvalue has real part Re A max ~ 38.5314476 for the Couette 
flow at Re = 5000 and fi = 40. 



in agreement with the general observation at the end of 
§11 A. Note that there is nothing wrong with unbounded 
growth rates for infinite values of /Lt, which is a common 
feature in many fundamental hydrodynamic instabilities 
in the short-wave limit, e.g. Rayleigh- Taylor instability. 
Of course, in the nonlinear setting one does not observe 
an infinite growth rates, which is suppressed by the non- 
linear effects, which can also be dissipative as will be 
shown in §111. Also with increasing /j, the eigenfunction 
structure becomes more complicated as we predicted in 
the above analytical study, e.g. for Re = 5000 and \i = 40 
the eigenfunction is shown in figure 4. More complete 



study of the stability picture for the Couette flow is out 
of the scope of this work and will be presented elsewhere. 



III. NONLINEAR TERMS EFFECTS 

As pointed out in the introduction, the usual argument 
made is that the linear terms in the NSEs are the only 
ones which produce energy and thus the nonlinear terms 
are energy-conserving and only mix energy between dif- 
ferent modes [4, 8, 9]. The claim made in [9] in the con- 
text of discussing shear flows with base state velocity field 
U is based on the consideration of the evolution of the 
kinetic energy E(t) = ||u|| 2 /2 of the disturbance velocity 
field u either over domains with periodic (i.e., compact 
domains) or homogeneous (zero) boundary conditions for 
disturbance, which leads to the classical Reynolds-Orr 
equation 



dE 
~dt 



Hivuir 



(Di 



(12) 



where Dij is the symmetric part of the base state ve- 
locity gradient tensor VU [50]. The linearity of (12), 
i.e. absence of the nonlinear terms effects, motivated the 
"transient-growth" idea of focusing only on the linear 
mechanisms. 

The problem with the above argument is that the 
"troublesome" shear flows are defined not on either over 
domains with periodic (i.e., compact domains) or homo- 
geneous (zero) boundary conditions for disturbance, but 
on domains with at least one semi- or unbounded dimen- 
sion and the boundary condition of boundedness of the 
disturbance field at infinity. Thus, let us re-derive equa- 
tion (12) taking into account that the boundary condi- 
tions on dn could be moved away to infinity. This leads 
to the following equation [51] 

dE 

~ ~dt z 



riiUipds + i 1 / rijUiU^j ds + v\\ Vu|| (13) 



on 
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on 



UjUjUiUi as 



+ t: / n^UjUiUi dx, 
2 A) 



where n is the normal outward (w.r.t. d£l) vector and 
part or the whole of dtt may be at infinity. If the do- 
main Q is unbounded, as in many applications, then the 
evaluation of the boundary terms in (13) becomes non- 
trivial, since it depends on the rate of the solutions decay 
in unbounded spatial directions [52], and therefore the 
nonlinear terms (cubic term in (13)) do not disappear, 
in general. In fact, there are no reasons to expect that 
these terms vanish in the Couette or other channel and 
pipe flows, since the disturbance, if it leads to an insta- 
bility, does not decay as it propagates from the entrance 
to infinity. Moreover, if we consider the flow on a semi- 
infinite domain, as discussed in §11, then there is a non- 
zero contribution of the nonlinear terms at this portion 
of dtt as well. As for open shear flows [53] the nonlinear 
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terms do not vanish since it is well-known that in order 
to construct a full set of eigenfunctions from the Orr- 
Sommerfeld equation, it is necessary to relax the homo- 
geneous boundary conditions at infinity to the condition 
of boundedness at infinity, |u| < oo, [10, 11, 12], be- 
cause the homogeneous boundary condition leads only to 
a finite set of a point discrete spectrum and therefore the 
inclusion of the continuous spectrum component is nec- 
essary for completeness. The completeness is of course 
important if one wants to study the question of stability 
rigorously. Therefore, one has to admit that the role of 
the nonlinear terms in the NSEs for shear flow problems 
is not a simple mixing of energy, but the disturbance en- 
ergy can be denerated/dissipated by the nonlinearities. 

Finally, another very important point is that the 
usual assertion that at early times the nonlinear effects 
are not important is not necessarily true since, for ex- 
ample, the linear term (Dij,UiUj) and nonlinear one 
~ \ Ion n j u j u i u i ds in equation (13) may become compa- 
rable even at small times either thanks to the very small 
values of Dij and/or the large cumulative effect of non- 
linear terms when integrated over dfl (possibly infinitely 
extended). In any case, no one has ever proved that this 
cannot occur for the flows in question! 

It should be kept in mind that the energy is a nonlocal 
measure of the fluid motion, while in reality we are inter- 
ested in the pointwise description, since we do not know 
much about the singular structure of the NSEs solution. 



IV. OTHER ISSUES 

For the purpose of the subsequent discussion discus- 
sion, it is convenient to treat the NSEs (1) as an infinite- 
dimensional ODE in the operator form in some Banach 
space X: 

^=Au + N(u), (14) 

where A is the linear operator, usually stationary, and 
N is the nonlinear operator. Historically, this approach 
allowed one to apply the dynamical systems methods of 
Dalekii & Krein [13], Yudovich [14] for nonlinear stability 
of the equilibrium solutions of (14). Equation (14) can be 
obtained from (1) by utilizing the Hclmholtz decompo- 
sition [15], i.e. through decomposing the solution vector 
field into the sum of its divergence- and curl-free compo- 
nents, e.g. L P {Q) = S p (fl) ® G p (fl), where S p (fl) is the 
Banach space obtained by closing the set of solenoidal 
vectors and G p (tt) is the Banach space obtained by clos- 
ing the set of gradients in the norm of L p (fl). Then 
one can introduce the Leray projector P, which projects 
any vector in L p (il) onto 5 p (fi), and thus yields (14). 
In particular, if the stability of a nontrivial stationary 
(time-independent) base state U(x) is studied, then the 



linear and nonlinear operators in (14) take the form 

Au = P [-U ■ Vu - u ■ VU + Re^Au] , (15a) 
N(u) = -P [V ■ (u <g> u)] . (15b) 

As follows from (15a), one can see the linear NSEs oper- 
ator A as a non-self-adjoint perturbation of the Laplace 
operator. 



A. Interrelation of stability notions 

The notions of stability or instability, as physically 
observable phenomena, were given by Lyapunov [16]. 
Namely, the equilibrium solution (base state) u e G X 
of (14) is said to be Lyapunov stable (sometimes called 
nonlincarly stable since (14) is nonlinear) if for any e > 
there exists a S > so that the initial conditions u G X 
and |juo — u e (to)||x < S imply that (i) there exists a so- 
lution u(i, £oj u o) and (ii) ||u(i,io,uo) — u e (i)||x < e for 
all t. Then, the notion of Lyapunov instability is simply 
the negation of the above definition of stability. While 
this negation formally does not require existence to hold, 
i.e. the solution seizing to exist is a particular form of 
instability, from the physical point of view one does need 
existence in X to get a sensible instability result [17]. 
Note that the linear stability is the Lyapunov stability of 
the linearized version of (14), u t = ^4u in the sense of 
the above definition. This should also be distinguished 
from the spectral stability, i.e. a formal notion obtained 
from the spectral problem, Av = Av, associated with the 
linearized version of (14), when the spectrum is in the 
left half-plane. 

There is a common view that linear eigenvalue analysis 
implies nonlinear Lyapunov stability, e.g. [4] who, in the 
context of their discussion of the transition to turbulence, 
refer to Sattinger [18] in order to claim that a connection 
between linear and nonlinear stability has been estab- 
lished. However, in that work Sattinger demonstrated 
stability using the Reynolds-Orr equation, that is in a 
non-pointwise norm, for compact domains only, while 
the shear flows are usually considered on domains with 
at least one extended dimension. In general, as one can 
gather from most of the fluid dynamics literature, it is 
used as a rule: if the spectrum of the linearized operator 
(15a) is in the left half plane then one has stability and 
the instability takes place if there are eigenvalues in the 
right half plane. Below we first discuss the interrelation 
of all these notions in finite dimensions, and then address 
infinite dimensions in §IVB. 

Example 1. First, consider a conservative system with 
the Hamiltonian H(p,q) = p 2 /2 + V{q), where the po- 
tential is quartic V(q) = q 4 /^. As one can immedi- 
ately see from the Hamilton's equations, q t = p and 
p t = —q 3 , the linear and nonlinear stability of its equi- 
librium, q = p = 0, definitions do not imply each other: 
this Hamiltonian system demonstrates nonlinear stability 
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because the energy function is concave, but its lineariza- 
tion around the origin, p = & q = p, produces a solu- 
tion growing linearly in time, i.e. it is linearly unstable. 
However, this example is spectrally stable; thus, spectral 
stability does not imply even linear stability, however the 
converse is true. 

One might think that the above is a peculiar feature 
of the Hamiltonian dynamics, but in fact one observes 
the same for dissipativc systems, e.g. for q% = P and 
Pt = -q - P ■ 

Another example, discussed in [19] (see also [20]) 
proves that linear stability does not imply nonlinear sta- 
bility: 

Example 2. Given the Hamiltonian 

H = \(p\ + «i) - (Pi + ?l) + \p-2(p\ - q\) - <7i<?2Pi- 

from the corresponding Hamilton's equations, we find that 
the origin, q = and p = 0, is linearly stable, but a one- 
parameter, t, family of solutions for this system 



In 



sin(i — t) 



Pi = \ 

qi = — \/2cos 



t-T ' 
(t-r) 

t-T 



Pi 
<12 



sin2(t - t) 
t-T ' 
cos2(t — r) 

t - T 



demonstrates nonlinear instability by blowing up at some 
finite time, t = r. 

Therefore, from the above examples we conclude that 

• linear instability does not imply nonlinear instabil- 
ity; 

• linear stability does not imply nonlinear stability; 

• spectral stability does not imply both linear and 
nonlinear stability, 

However, spectral instability implies nonlinear instabil- 
ity in finite dimensions based on the classical Lyapunov 
indirect method [16]. 



B. Infinite dimensionality effects 

In infinite dimensions, the situation is even richer, be- 
cause the latter property - the Lyapunov indirect method 
- does not hold in general. As mentioned in §IVA, the 
usual approach is to base the conclusion of stability or 
instability of (14) on the spectrum a(A) of the linear op- 
erator A: if all eigenvalues lie strictly (i.e. no eigenvalues 
are on the imaginary axis) in the left half of the complex 
plane C, then the zero solution is stable, while if there 
are eigenvalues in the left half-plane of C, then the zero 
solution is unstable. 



Let us discuss the question of linearized stability sys- 
tematically [54] . The solution of the linearized version of 
(14) which can be written symbolically as 



u(t) = u(< ) e 



A(t-t ) 



(16) 



For a Banach space X of finite dimension it is well known 
that if sup A G <j(A) is < 0, then ||T(i)|| decays exponen- 
Re A 

tially [21]. This behavior is a consequence of the fact 
that linear operators in finite-dimensional Banach spaces 
have only point spectrum. Since this is not the case in 
infinite-dimensional Banach spaces one does not expect 
this result to be true in general. From a formal point 
of view, the exponent e At makes mathematical sense if 
the operator A is bounded, i.e. ||At||x < M||x||x (in 
future we drop the index denoting the space X in which 
the norm applies unless it becomes important to distin- 
guish), so that 



e At = Y J A k t k /k\, 



(17) 



k=0 



is well-defined, and from the spectral mapping theorem 
[22] one can connect the spectra a(A) and a(e At ) via 



cr(e 



AL\ 



(18) 



The convergence of (17) allows one to get a rough esti- 



mate, He^H < J2 \\A\\ k t k /k\ 

k=0 

the knowledge of the spectrum a(A) allows one to get 
sharper estimates: in particular, if ||e || < q, the spec- 
trum <7 (A) lies in the interior of Re A < log q. If, however, 
the operator A is unbounded, which is usually the case in 
hydrodynamics, then e At makes only symbolic sense and 
the equality is not true in general [23] (see also Lax [24], 
pp. 434-439) and, in fact, the spectral mapping property 
(18) does not hold, i.e. 



\\A\\t 



t > 0. However, 



a(e At ) D 



(19) 



from where it follows that one cannot conclude on stabil- 
ity of (14) based on the knowledge of spectrum of cr(A). 
For unbounded operator the property (18) takes place 
only if it is an infinitesimal generator of an analytic semi- 
group, the convenient test for which is via proving that 
the operator is sectorial (cf. Henry [25]), i.e. if its eigen- 
values are contained in a cone sector with an apex angle 
< 7r in the left-half of the complex plane. 

These known facts can be most conveniently formu- 
lated with the help of two notions: spectral bound, 
r{A) = sup {Re A| A € o~(A)}, and the growth abscissa 



,At\ 



Then, 



of the semigroup, u>{A) = lim t 1 log (| 

t— »oo 

in finite dimensions it is true that u>(A) = r(A), which 
is implicit in all studies of hydrodynamic stability. In 
the infinite-dimensional case, this equality is not true in 
general, as will be discussed below. However, some op- 
erators in certain spaces do possess this property, such 
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as the Stokes operator is known [14] to generate an an- 
alytic semigroup in any L p space. Thus, one can expect 
that if (15a) is not a strong enough perturbation of the 
elliptic operator, which also depends on the boundary 
conditions, then the operator should still be a generator 
of an analytic semigroup. 

On the other hand, there are many stability problems 
where the diffusive mechanisms, as controlled by elliptic 
operators (the Stokes, biharmonic, etc.), may be negli- 
gible, and thus u>(A) ^ r(A). For example, hyperbolic 
PDEs exhibit such behavior, as can be illustrated with 
the following simple example due to Renardy [26] 

Example 3. Consider the following hyperbolic equation: 

u u = u xx + u.yy + e iy u x , (20) 

with 2tt -periodic boundary conditions in both directions. 
Introducing v = Ut, the semigroup operator A associated 
with (20) is A(u,v) = (y, u xx + u yv + e w u x ), which in 
fact generates a Co-semigroup in the space H 1 xL 2 for the 
pair (u,ut). The analysis [26] demonstrates that r(A) = 
0, while J(A) = 1/2. 

In general, it is well-known after the work of Zabczyk 
[27] that for any two real numbers ui s < u>q there exists 
a strongly continuous semigroup {T t } t > on a Hilbcrt 
space, such that ui s = sup {Re A; A £ o-(A)} and |T t | = 
e w °*, t > 0, where a (A) denotes the spectrum of the 
generator A of the semigroup {T t }t>o- 

The above most systematically can be understood 
in terms of classification of semi-groups [21], namely 
parabolic versus hyperbolic, as resulted from the appli- 
cation of the abstract semi-group theory to evolution 
partial differential equations. In the parabolic case the 
operators are the infinitesimal generators of an analytic 
semi-group, while in the hyperbolic case one can expect 
the behavior as in the above example 3 due to Renardy 
[26] . While behavior of example 3 is common for Hamil- 
tonian systems, in which the eigenvalues are distributed 
symmetrically w.r.t both axes in the complex plane and 
the placement of all eigenvalues on the imaginary axis is 
a necessary but not a sufficient condition for stability of 
an equilibrium, equation (20) is clearly non-conservative, 
which suggests that such pathologies in the linear stabil- 
ity picture are common for dissipative systems too. 

First of all, the existence and (in)stability of the so- 
lutions of the full nonlinear system (14) are most con- 
veniently (at the current level of mathematics) studied 
by transforming it into the integral form via Duhamel's 
formula 

u(t) = T(u) = u{0)e At + J e A ^- T) N(u(t)) dr, (21) 

o 

which is not only used to solve (14) using fixed point the- 
orems for the mapping 7~(u), but also clearly illustrates 
the importance of both the linear terms - the semigroup 



{e At }t>o ~ an d the nonlinear terms in infinite dimen- 
sions. Indeed, on any relevant to this discussion Banach 
space X the nonlinear terms will be unbounded in view of 
the presence of unbounded operator of differentiation in 
(15b), which can be compensated only by the smoothing 
effect of the semigroup e A ^~ T \ The latter can happen, 
of course, only if the semigroup is "nice" enough. This 
illustrates the importance of consideration of both lin- 
ear and nonlinear effects when studying the questions of 
stability. Formula (21) is the basis of all rigorous sta- 
bility studies, such as due to Krein [13], Henry [25], and 
Yudovich [14]. 

Thus, it is not surprising that infinite dimensionality 
may lead to pathologies, which are not possible in fi- 
nite dimensions. For example, in the context of elastic- 
ity problems it is known that the positive definite sec- 
ond variation of the energy (and thus the correspond- 
ing system should be nonlincarly stable according to the 
Dirichlet theorem, valid in finite dimensions only) does 
not guarantee the stability in view of the possible pres- 
ence of an infinite number of unstable directions [28] . 

Another complication arising due to infinite dimen- 
sionality is the norm-dependence of stability criteria - the 
issue which has been understood for a long time [14, 29]. 
In the finite-dimensional case this difficulty cannot arise 
since all norms in finite-dimensional Banach spaces are 
equivalent. 

Example 4. One of the simplest examples of this sub- 
tlety is due to [H] and represents a linear PDE: 

du du 
dt dx' 
u(0, x) = 4>{x), 

the unique solution of which is u(x,t) = <p (x e*). Thus 
one can express the L p '-norm of the solution derivatives 

via ^tt± =e {k ~ p ~ 1)t U (k) \\ TB ,™, and therefore 

one has (a) asymptotic stability in L P (R) for 1 < p < 
oo, (b) Lyapunov stability mL°°(R), and (c) exponential 
instability in Sobolev spaces W k ' p (M.) with k > 1, p > 1 
or k = 1, p > 1. 

Finally, quite often by "nonlinear stability" one under- 
stands "energy stability" , which is not quite correct, since 
nonlinear stability is a stability in the Lyapunov point- 
wise sense, while the energy norm is a global measure 
[55] . In order to make a pointwise sense of the energy- like 
norms, so that the function belongs to the Holder space 
of fc-times continuously differentiable functions f E C k , 
one needs energy-like bounds for the solution derivatives, 
that is / should belong to the Hilbert space H k with 
high enough index, so that one can employ the Sobolev 
inequality; e.g. in the two-space dimensions: 

\f\ ck <C\\f\\ Hk+Sl s>l. 

From here it follows that establishing bounds on the 
usual kinetic energy norm H° is not enough to assert the 
bounds on the solution, and thus to claim the stability 
of the solutions in the Lyapunov sense. 
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C. Non-normality and covariance 

Starting with Trefethen ct al. [3] the following two- 
dimensional model is very popular [4] for illustrating 
what is "supposed" to happen in the transition problem 



du 



Au + F(u), A 



~CL\ 1 

-a 2 



(22) 



where A is a non-normal operator in a sense that AA*^ 
A* A, where A* is the adjoint operator, with ai,2 ~ Re -1 
being small exponents determining the time-evolution, 
and F(u) is a nonlinear energy conserving operator and. 
Since the linear operator A is non-normal, then its eigen- 
vectors corresponding to the eigenvalues Ai = —a\ and 
A2 = — a.2 are almost parallel for Re 3> 1, which in 
the case used by Trefethen et al. [3], a± = Re -1 and 
a.2 = 2 Re^ 1 , are: 



(23) 



In the limit Re — > 00 the matrix A becomes the non- 
trivial Jordan block, when the solution u grows alge- 
braically. Due to this non-normality for finite Re, the 
system u = Au experiences a transient growth for t = 
O(Re), e.g. for the initial conditions u(0) = [0, 1] T : 





' 1 " 




1 


ei = 





, e 2 = 


-Re' 1 



u(t) 



Re 






' -Re ' 


2t 


e «<= + 


1 


e 



t + O 



1 



— + o 

Re ^ ^ 



t 

R, - 



(24) 



B = 



Next, the key dynamic feature of all these models is 
the presence of a finite-amplitude instability; in the case 
of Trefethen et al. [3] the nonlincarity is ||u||Bu with 

j ^J, which has three stable nodes (one at the 

origin) and two saddle points, while Manneville [30] used 

F(u) = (uv, —u 2 ) T , which has two stable nodes (one is at 
the origin) and a saddle. Manneville used example (22) 
to advocate that the non-normality of the linear opera- 
tor is not in itself responsible for the by-pass transition, 
while Trefethen et al. [3] insisted that the transition is 
"essentially linear" . While the form of Trefethen's non- 
linearity was criticized [31] as not relevant to the NSEs, 
it is clear that it cannot be obtained as a result of some 
kind of reduction from the NSEs, since the square root 
is not Taylor expansion-friendly. 

If one expects that transition is a fundamental physi- 
cal phenomenon, then it must be accounted in a covari- 
ant (coordinate-free) manner, that is its understanding 
should be independent of the coordinate system. Thus, 
let see if model (22) possesses this property. 

In this subsection we consider the linear part of (22) 
recalling some elementary facts, while in the next section 
we discuss the fully nonlinear dynamics. Even though 



the original matrix A in (22) is non-normal, there ex- 
ists a transformation which makes this linear operator 
in (22) normal. Clearly, normal matrices are a subset of 
the diagonizable matrices and the relation between nor- 
mal and diagonizable matrices is well-understood [32]. 
The notion of diagonizability is an intrinsic notion (that 
is, independent of a coordinate system) as opposed to 
to non-normality, since the latter notion depends upon 
a system of coordinates. In terms of genericity notion, 
cf. Wiggins [33], i.e., informally, how common are the 
physical systems with non- normal operators, it is obvi- 
ous that non-normal matrices form a dense subset in the 
2 

set C n of all n x n matrices. Further, we will comment 
on the infinite-dimensional operators, but first consider 
the finite dimensional case as relevant to the discussion of 
model (22). In fact, this clear understanding of which no- 
tions are intrinsic or not shades some light on the claims 
that the non-normal linear operators is the key to transi- 
tion. Appealing to the standard matrix theory [34] , let us 
consider a linear transformation A on the Euclidean space 
E™, which maps the vector x into x': x' = Ax. If matrix 
T consists of the basis vectors t 1 , . . . ,t™, then x = Ty 
and x' = Ty' , where y and y' arc components of the cor- 
responding vectors x and x'. Note that det T^O since it 
forms a basis. Then, the relation between y and y' is sim- 
ply y' = T" 1 A T y. Matrices A and T^ 1 AT are similar, 
and T~ J AT can be made diagonal (and thus normal) if 
and only if A has n linearly independent eigenvectors [56] 
(in particular, if A has distinct eigenvalues, but not nec- 
essarily). Of course, if the matrix A is originally normal, 
then it will have orthogonal eigenvectors and thus T will 
be orthogonal, which means that normality will be 'pre- 
served'. Anyway, non- normal matrices may be reduced 
to normal ones under appropriate change of coordinates. 
Since we looking for a covariant description of physical 
phenomena, i.e. independence of our understanding of 
the phenomena of a coordinate system, then it becomes 
clear that non-normality has nothing to do with the co- 
variant understanding of the fundamental cause of the 
transition. While the transient growth effects may be 
important exactly as "transient" effects, which can be 
most clearly seen through the singular value decomposi- 
tion, on the time scales greater than the transient growth 
time, e.g. t > Re" 1 in (24), their effect is not relevant in 
any coordinate system. Becides these general remarks, 
one should also point out from the positions of control 
theory that not all non-normal operators can lead to sig- 
nificant transient effects, but only those which possess 
some high sensitivity subspaces in their domain of defi- 
nition together with the presence of disturbances in this 
subspace. 

The above discussion for operators in finite- 
dimensional spaces can be translated to the case of 
infinite-dimensional operators, as we discuss here in the 
context of PDEs. Clearly, the linear operator in the 
NSEs (15a) is non-normal for non-trivial base states 
U(x). The degree of non-normality of (15a) increases 
as the Reynolds number increases since the self-adjoint 
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part of (15a) [57], that is the Laplacian, becomes less 
dominant. Most clearly, this probably can be seen in a 
spectral space after the Galerkin projection. In any case, 
the non-normality of the partial differential operator 
leads to a transient in time growth in the same way as a 
non-normal matrix-operator, which is a self-evident fact, 
but nevertheless quite a number of works were devoted 
to illustrate this behavior in the context of the NSEs 
[5, 35]. 



D. Finite-amplitude effects 

From the nonlinear covariant analysis viewpoint, it is 
tempting to perform the normal form analysis of (22) in 
order to reveal the universal behavior of (22). However, 
the local nature of the normal form analysis, i.e. inability 
to get the global bifurcation picture, does not allow one 
to capture the finite-amplitude instabilities inherent in 
(22): e.g. in the non-resonant case of a\ and a-i in (22), 
the normal form is simply 



dx 
~dt 



Jx, J 



-en 

-02 



and in the resonant case, ai = 2a\ is 





dx 
di 



Jx 



ax\ 



, a G 



Both the above normal forms do not exhibit finite- 
amplitude instabilities. This fundamental inability of the 
normal form analysis to capture the global bifurcation 
picture can be understood with the following simple ex- 
ample, which also leads to some further insights. 



Example 5. Let us consider the following equation 



dx 



-ex + x' 



(25) 



which describes transcritical bifurcation and exhibits a 
finite- amplitude instability with the critical amplitude 
x c = e. With the goal to reduce this equation to a linear 
equation 



dy 

di = ^ 



(26) 



let us introduce a transformation, x = y + F(y), where 
F(y) is a function to be determined. Substitution of this 
transformation into (25) gives 

^ [1 + F'(y)] = -e[y + F(y)} + [y + F(y)} 2 . (27) 

In the standard treatment of the local normal form anal- 
ysis (e.g. [33]), one assumes that \F'(y)\ <C 1 and thus 
inverts [1 + F'(y)] approximately, i.e. [1 + F^y)]^ 1 ~ 
1 — F'(y). In fact, one does not need this approximate 



procedure in order to get an equation for F(y): instead, 
one multiplies (26) by F'(y) and subtracts from (27): 



dy 
dt 



-ey - eF(y) + eyF'(y) + [y + F(y)} 2 , (28) 



that is, F{y) is determined from 

e[F-yF'} = [y + F} 2 , 



with the result 



F(y) 



y 



y, 



(29) 



(30) 



where a / for the transformation x = y + F(y) to be 
non-singular. Thus, for x ^ e the transformation x = 
y + F(y) "straightens out" the trajectories of the original 
system (25) and produces a trivial dynamics as in (26), 
but is singular for x = e. 

The latter property, i.e. singularity of the transfor- 
mation which linearizes the dynamics globally except for 
at the threshold amplitudes can be used to identify the 
presence and location of the finite-amplitude instabilities 
in the appropriate configuration (phase) spaces. 

The indicated in the introduction possibility of finite- 
amplitude instability nature of the transition suggested 
a search for finite-amplitude solutions in the shear flows. 
For example, Nagata [36] found a finite-magnitude peri- 
odic solutions in the Couette flow, which coexists with 
the linear Couette profile, appealing to the concept bi- 
furcation from infinity [37] . The travelling wave-like so- 
lutions are known both in the context of the pipe [38] 
and plane Couette [39] flows. However, their "strongly 
unstable character", since they are all saddle points in 
phase space [40], is also well-known [41]. The latter work 
also advocates that the turbulent state in pipe flow cor- 
responds to a chaotic saddle (unstable aperiodic orbit) 
in state space. The idea is that travelling wave solutions 
presumably constitute a 'skeleton' about which compli- 
cated time-dependent orbits may drape themselves tem- 
porarily before falling back to the laminar state [40] . 

Waleffe [31] advocated the idea of exploring the size 
of the domain of attraction of the finite-amplitude solu- 
tions following the original thoughts of Orr and Thom- 
son. Also, Trefethen et al. [3] conjectured that the non- 
normality of the linear operator shrinks the size of the 
attraction basin and, in fact, the threshold amplitude 
scales as i?e -7 . 

The above dynamical systems approach is based on the 
finite-dimensional view of the NSEs dynamics. The lat- 
ter is usually justified [40] by the argument that the mo- 
tion of a viscous fluid in a finite domain is always finite- 
dimensional which is due to the viscous cutoff of fine 
scales with reference to [42]. However, the troublesome 
flows always have at least one unbounded dimension, 
which undermines the above logic and poses the question 
on the possible crucial effect of infinite-dimensionality. 
Next, as currently understood, all these finite-amplitude 
solutions proved to be unstable. Moreover, it is very 
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likely that if the transition is indeed a finite-amplitude 
instability, the attracting solutions in fact are intrinsi- 
cally time-dependent or even chaotic and occupy some 
subset of the phase space. Thus, there is little hope to 
find the corresponding solutions analytically. 

V. DISCUSSION 

The purpose of this note was to bring a number of 
important theoretical issues to the attention of the flu- 
ids community, solving of which may help to make some 
progress on a reasonable qualitative understanding and 
quantitative prediction of the experimental observations, 
which are still lacking. In particular, as follows from §11 
and §111, exploring the effects of the domain type and of 
the energy conservation by nonlinear terms are the first 
natural issues to address systematically. 

The authors also see two other possible ways of explor- 
ing the problem of transition rigorously. The first way is 
to develop a semi-local nonlinear stability theory, which 
could be a natural generalization of the stability theory 
due to Krein, Yudovich, and Henry by inclusion of the 
finite-amplitude instability picture in consideration. The 
second possible avenue is to gain insights into the geomet- 
ric structure of the equations, which allow one to identify 
the possible regions of attraction in the phase space and 
thus to establish the stability picture in large (i.e. both 
local, global, finite-amplitude, etc.). Namely, one first 
explores the phase space structure of the Hamiltonian 
approximation (ideal fluid) and, if this is not sufficient 
for explaining the instability picture, one adds dissipa- 
tive effects (viscosity). In particular, this would allow 
one to understand if the fundamental basis of the tran- 
sition phenomena is Hamiltonian or intrinsically due to 
dissipative (viscous) effects; see also [17]. The reader can 
easily appreciate the effectiveness of this approach with 
the help of, say, the Takcns-Bogdanov system [43] : 

x = n-x 2 = -V x , V{x) = -fix + x 3 /3. (31) 



Namely, locating maxima and minima of the potential 
V(x) (or alternatively of the Hamiltonian) indicates the 
stability picture. This idea, of course, is very old and goes 
back to Lagrange, Dirichlet and other classics. Next, 
one adds dissipative effects [44] and observes how the 
stability picture changes. The infinite-dimensional case 
is, of course, more complicated, but some progress has 
been done in this direction as well [17]. 

Finally, while all the previous fluid mechanics experi- 
ence over the last few centuries suggests that the NSEs 
arc the adequate description of fluid motion, strictly 
speaking one cannot discard the chance that the NSEs 
do not describe the subtle transition to turbulence phe- 
nomena. One can name many reasons for which the 
NSEs may turn out to be inappropriate for modeling 
of the transition. For example, since the transition is 
a phenomenon presumably very sensitive to initial and 
boundary conditions, then it should be very sensitive to 
the details of the equations. Since Newtonian fluid de- 
scription is an approximate one, then in the conventional 
NSEs for incompressible fluid we discard by hyperbolic 
effects common for non-Newtonian fluids. Same can be 
said about the incompressibility approximation. While 
one might argue that these effects may have an influence 
only at some marginal time and spatial scales, the his- 
tory of hydrodynamics knows a number of fundamental 
examples when small effects affect the flow in the large 
(Prandtl's boundary layer theory, for instance). 
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